Gravity wave analogs of black holes 
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It is demonstrated that gravity waves of a flowing fluid in a shallow basin can be used to simulate 
phenomena around black holes in the laboratory. Since the speed of the gravity waves as well as 
their high-wavenumber dispersion (subluminal vs. superluminal) can be adjusted easily by varying 
the height of the fluid (and its surface tension) this scenario has certain advantages over the sonic 
and dielectric black hole analogs, for example, although its use in testing quantum effects is dubious. 
It can be used to investigate the various classical instabilities associated with black (and white) holes 
experimentally, including positive and negative norm mode mixing at horizons. 
PACS: 04.70.-s, 47.90.+a, 92.60.Dj, 04.80.-y. 



I. INTRODUCTION 

One of the most fascinating predictions of Einstein's 
theory of general relativity is the potential existence of 
black holes - i.e. space-time regions from which nothing 
is able to escape. Perhaps no less interesting are their 
antonyms: white holes (nothing can penetrate). Both 
are described by solutions of the Einstein equations and 
are related to each other via time-inversion, see e.g. 

As it is well-known, these objects feature many novel 
properties: For example, for orbits sufficiently close to 
the horizon (i.e. for r < 3Af ) one observes 0] an inversion 
of the centrifugal acceleration. 

Rotating black holes as described by the Kerr metric 
admit unstable modes under certain conditions, i.e. solu- 
tions of the wave equation growing in time without any 
bound, see e.g. pLpJ- This phenomenon is related to the 
mechanism of superradiance || which allows one to ex- 
tract energy from the rotation of the Kerr black hole, cf. 

White holes are unstable || to exponential build-up of 
energy on the white hole Cauchy horizon on the classical 
level, as well as on the quantum level [Q-f)). 

The presence of both Cauchy and particle horizons 
(white and black hole horizons), such as in the interior 
of a Reissner-Nordstrom metric, can have further insta- 
bilities, see e.g. JTo| . 

Another striking effect is the evaporation |ll|] of black 
holes due to quantum effects. This observation can be 
interpreted as a confirmation of their thermodynamical 
interpretation [ fl2[ relating purely geometrical quantities, 
such as surface gravity and surface area, to thermal prop- 
erties, such as temperature and entropy. 

Fortunately it seems unlikely that one can observe 
black holes in the laboratory (see, however, e.g. 
Analogs, which obey similar equations of motion to fields 
around a black hole raise the possibility of demonstrat- 
ing some of the most unusual properties of black holes in 
the laboratory. This is the basic idea of the black and 
white hole analogs (Dumb holes) originally proposed by 
Unruh in Ref. [B. The sonic analogs established there 
are based on the observation that sound waves in flowing 



fluids are (under appropriate conditions) governed by the 
same wave equation as a scalar field in a curved space- 
time. The acoustic horizon, which occurs if the velocity 
of the fluid exceeds the speed of sound within the liquid, 
acts on sound waves exactly as a black hole horizon does 
on, for example, scalar waves. 

After the original proposal in Ref. H] the sonic analogs 
have been the subject of several investigations, see e.g. 
Jl5| [TtJ . Although the kinematics of the waves propagat- 
ing within the black and white hole analogs are governed 
by the same equation as those in a curved space-time, the 
dynamics of the effective metric itself are not described 
by the same laws as gravity (i.e. the Einstein equations) 
in general^]. 

In this way the analogs allow one to separate the dy- 
namical effects of gravity (following from the Einstein 
equations) from more general (kinematic) phenomena, cf. 
©[§■ 

In addition to the sonic analogs there exist proposals 
for black or white hole analogs based on the propagation 
of light in dielectric media (instead of sound), see e.g. 
p0| p4| , and of other sorts of waves in for example liquid 
Helium 3, see e.g. These scenarios avoid some of the 
difficulties associated to the sonic analogs but can have 
other problems. 

The challenge in making such analogs to horizons is 
in preparing a medium in which the waves are stopped 
from propagating out from some region. In the analogs 
where the flow of a medium is used to drag the waves 
at a velocity corresponding to the velocity of the waves, 
one requires a sufficiently low velocity that the experi- 
ment could be contemplated. The speed of sound de- 
pends on the equation of state p = p(g) only and there- 
fore is hard to adjust by external parameters. In the case 
of the analogs based on light, the velocity of the (quasi) 
photons is determined by the effective permittivity and 
permeability of the medium, which can also be hard to 



* There are, however, possibilities to reproduce the Einstein 
equations even in non-gravitational systems, see (l8|. 
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manipulate, and especially hard to make a sufficiently 
low group and phase velocity of light (cf. 22 2^,|25|]). 

Consequently we were led to look for another kind of 
waves traveling at a velocity that can be controlled more 
easily. One promising candidate is gravity waves (sur- 
face waves) in a shallow basin filled with a liquid. As 
we shall see in the following, long gravity waves within a 
flowing fluid are also governed by the same wave equa- 
tion as a scalar field in a curved space-time. In addition, 
the speed of the long gravity waves can be adjusted very 
easily by varying the depth of the basin. Furthermore 
the fluid flow in such a basin is easily manipulated. Be- 
cause of the low velocity of these waves, quantum effects 
would not be observable, but many of the classical fea- 
tures of black holes (including the positive and negative 
norm mixing at the horizon which is closely related to 
the quantum evaporation effects) could be investigated. 
Furthermore, as we shall see, the dispersion relation of 
these waves at high wavenumbers can also be easily ma- 
nipulated, allowing easy investigations of the effects of 
such changes on horizon effects. 



II. THE MODEL 

We shall begin with the simplest form of the model, in 
which we assume a shallow liquid over a flat, horizontal 
bottom. Furthermore, the forces on the liquid will be 
assumed to be such that they allow for a purely horizontal 
stationary flow profile resulting in a constant height (i.e. 
horizontal surface) of the liquid. Later we shall relax both 
of these assumptions. In addition, we shall assume that 
the liquid is viscosity free, incompressible and irrotational 
in its flow. 




FIG. 1. Picture of a gravity wave in the basin and the 
relevant parameters. The assumed relation of the dimensions 
for long gravity waves 5h <C /ib <C A is not reproduced for the 
sake of conciseness. 

In such a case the density of the liquid remains con- 
stant (g = const) and in terms of its local velocity v the 
equation of continuity assumes the simple form 



V -v = 0. 



(1) 



If we neglect the viscosity of the fluid its dynamics are 
governed by the non-linear Euler equations, see e.g. pa] 



— = v + (v-V)v = -+g + L , 2 

at g g 

with p denoting the pressure and g = —g e z the gravita- 
tional acceleration; and / = — pV||V" is some horizon- 
tal force necessary to establish the stationary horizontal 
flow; cf. Fig. |. 

For an irrotational flow profile V x v = we may 
simplify the Euler equation (^) via (v ■ V)i? = V(v 2 )/2 
and introduce a velocity potential v = V$ arriving at 
the Bernoulli equation 



$+i(V$) 2 = --- ff Z-T/H 

2 g 



(3) 



The boundary conditions are that the vertical velocity 
must be zero at the bottom of the tank, the pressure zero 
at the displaced surface, and change in the height of the 
fluid determined by the vertical velocity (cf. Fig. [l]) 

v ± (z = 0)=0, v ± (z = h) = ^ = h+(vV)h, (4) 



and 



p(z = h) = 0. 



(5) 



Now, let us consider perturbations Sv to a background 
flow «b (which is assumed to be stationary, irrotational, 
and horizontal) corresponding to small vertical displace- 
ments 5h of the height of the fluid, h. The background 
flow will be assumed to obey 



i.e., tiB — const, and 



V|| -vb = 0, 



2 g 



(6) 



(7) 



where p B will be given by g(h B — z). 

We shall assume that the velocity perturbations are 
also irrotational, so that they are given by a potential, 
<5<I>. The perturbations of the Bernoulli equation are 
given by 



S$ + w| • V||<5$ = -— . 



(8) 



The boundary condition (|^) for the pressure together 
with pb = ff(^B — z) imply 

5p(z = h B ) = ggSh, (9) 

and similarly for the vertical velocity 

6v ± (z = h B ) = Sh+ (v B -V\^ Sh, (10) 

as well as 5v J -(z = 0) = 0. 

It is useful to expand the perturbation potential <5<f> 
into a Taylor series 
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5<$>(x, y,z)=^2 - r «y$ (n) (x, y) 



(11) 



The boundary condition in Eq. ([3]) implies <5$(i) = 0. 
Another constraint arises from the equation of continuity 

(ED 



|d<J> (0) + <5$ (2) = . 



(12) 



and so on for larger values of n. Wc assume that the 
wavelength A of the perturbation is much longer than 
the depth h B . In this case the higher-order terms in 
the Taylor expansion (^Tj) are suppressed by powers of 
/ib/A <C 1 since we have V 2 = 0(1/ A 2 ). Keeping only 
the two lowest (non-trivial) terms in Eq. ([jj]) we find that 

8v ± (z = h B ) = -h B VlS<P {Q) . (13) 

This enables us to combine Eqs. (||), (||), and ( ^0| ) 



d 

at + " 



B 



V I <5$ (0) - £//t B V|<y$ ( o) = . (14) 



This wave equation, however, equals the Klein-Fock- 
Gordon (KFG) equation 



□ 5$ 



1 



(0) 



-^y^ggTduS^)) =0, 



(15) 



with the effective metric (remember V11 • v B = and 
g/iB = const) 



where ij go over values 1 and 2 (Einstein summation con- 
vention), and represent the coordinates within the bot- 
tom of the container. The equation of continuity ([!]) is 
now 



(19) 



with r\ = det(rjij) and, assuming irrotational flow, the 
Bernoulli equation (||) becomes 

$ + r) ij (di$)(dj$) + (9 2 $) 2 = -~ Q - V{x\z) , (20) 



with 



and 



d z <S>. 



(21) 



(22) 



Here the potential V(x l ,z) already includes the gravita- 
tional acceleration - in contrast to the potential V" used 
in the previous Section. 

The surface of the liquid, defined by z — h(x^), is where 
the pressure goes to zero and obeys, cf. Eq. (K 



h + v'dih = v z 



(23) 



where v % and v z are evaluated at the surface. 

Let us now expand these expressions in powers of the 
vertical height z above the bottom at z = 0. The velocity 
potential <&, the metric Ty*- 7 , and the potential V can be 
written as 



flcff = 



gh B l 



(16) 



Except for replacing the velocity of the gravity waves 
V 'gh B by the speed of sound it is exactly the same ef- 
fective metric as for the sonic analogs. Calculating the 
inverse q°JI of the effective metric one obtains 



off 

1 00 
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'gh B 



(17) 



As one would expect, the condition of an ergosphere 
0oo = is fulfilled where the velocity of the fluid 
equals the speed of the gravity waves \/gh B . 



III. ARBITRARY BOTTOM AND HEIGHT 

Let us now relax the previous assumptions that the 
bottom and the background flow surface are both flat 
and parallel. We introduce arbitrary coordinates on the 
bottom and define a vertical coordinate z as orthogonal 
to the bottom of the container and geodesic. The spatial 
metric can always be cast into the form 



dr 1 = dz 2 



Vij 



dx l dx 3 



(18) 



*(*«,*) = $(„)(*«) + y<f(2)(z l ) +0(z 3 ) , 
^{x\z) = ^{x l ) + zrf^(x l ) + 0(z 2 ) , (24) 



and 



V{x\z) = V m (x l ) + zg z (x l )+0{z 2 ) 



(25) 



where we already have incorporated the boundary con- 
dition v z (z — 0) = and introduced the gravitational 
acceleration perpendicular to the bottom g z = Vm . Sim- 
ilarly we obtain for the pressure p 



p=(h-z)p (1) +0([h-z] 2 ). 



(26) 



In analogy to the previous Section we assume the height h 
of the fluid to be much smaller than the horizontal length 
scales on which the features of the flow profile (e.g. rfi , 
V , and p) change significantly - such as the wavelength 
A. In this long- wavelength limit the higher-order terms 
of the above Taylor expansions are supressed by powers 
of h/X and thus can be neglected. 

The continuity equation (fy) enforces again 



1 



(2) 



(2) 



(27) 
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Evaluated at z = h the equation fl23) for the height in 
terms of the velocity reads 



h + if* [d.^(x k , h)] djh = d z <S>(x k ,h) , 
which again to lowest order in z = h becomes 
h + ^(d^djh = - W|#( ) + C(h 2 ) , 
or, equivalently, 



0. 



(28) 
(29) 

(30) 



which can be interpreted as an effective conservation law 
h + V|| • {hvU) = 0. 

The Bernoulli equation transforms into 



*(o) + 2 G(o) + ZT $)) (ft*(o))(»i*(o)) 



(h-z)- V (0) - zg z + 0{h 2 ) 



(31) 



(32) 



From the terms linear in z we may infer 

Pj f=9z+\rt^{d^ {0) ){d^ {0) ), 

We can define an effective gravitational acceleration as 
1 



S = 0»+ O*#)( 5 i*(O))(0i*(O))- 



(33) 



Note that the eigenvalues k of rj 1 ^ with respect to rj^L 

(34) 



i7 11 

'I(l) x 3 — K 'l(0) x 3 ' 



are just twice the inverse of the principle radii of curva- 
ture of the surface over which the fluid is flowing. Thus 
the extra terms in Eq. ([53]) just represent the vertical 
centrifugal forces on the fluid travel over this curved sur- 
face, and in all of our following investigations will be 
negligible. 

Let us assume again that we have a stationary back- 
ground flow which obeys these equations, and we are 
interested in perturbations around this flow. The per- 
turbation equations then are (using v B = Wfn\dj$f \ at 
z = 0) 

1 



Sh + di (JW^v B Sh) 



and 



£$ ( o) + v B di6<f> {0) = -gSh. 
We can combine these equations to get 



(35) 



(36) 
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Wj)8j »(0)=0. (37) 



This again is a KFG equation with a metric given by 



- — 



1 



J B "B 



ghBV? 



(0) 



(38) 



where both g and /ib can depend on the coordinates 
x k . We can thus sculpt the effective metric within which 
these waves flow both by altering the velocity of the back- 
ground flow, by changing g from place to place (primar- 
ily by sloping the bottom of the tank), or by altering 
the height of the background flow. (Of course the back- 
ground /ib is determined by the background flow, g and 
the potential V.) If one has only the gravitational field as 
a force on the fluid, the slope of the bottom can be used 
to generate a potential V, and also, with more severe 
slopes, to change the value of g from place to place. 

The ergoregion is defined as the zone where, in order 
to be travelling at less than the velocity of the wave in 
the rest frame of the fluid, one cannot be standing still 
in the lab frame. This is the region where the velocity 
of the fluid is higher than the local velocity of the wave, 
VWb ■ 

If we assume that the bottom is flat (so that rfi = 5^ ) , 
and that the flow is driven by changes in /ib, we have 



-v B + gh,B = const . 



(39) 



But \ZgliB is the speed of the gravity waves. Far from 
the ergosphere, the velocity vb is small, so the we have 
const = ghoo. Thus at the erogosphere, where the veloc- 
ity of the fluid is the velocity of the waves, we have 

h = ^h 00 . (40) 

Similarly, if we assume that we have a sloping bottom de- 
signed so that the fluid maintains a constant height /is 
throughout, we would again obtain that the ergosphere 
should be at a point such that V{x l ) — = gh B /2. If 
this potential arises purely from the gravitational poten- 
tial due to the slope of the bottom, we must have that the 
bottom would have to be at a height of h/2 lower than 
at infinity. I.e., it does not take much of slope to the 
bottom to create the conditions necessary for an analog 
black hole ergosphere to form. 

We note that in the above, we have only kept terms to 
lowest order in z or /i. The validity of this approximation 
is essentially that all horizontal derivatives have scales 
which are much larger than h, the height of the fluid. 
Furthermore, the condition that we need only retain the 
lowest order in the metric rfi is also that the curvature 
of the bottom of the tank be on scales which are long 
with respect to h. I.e., we are in a "shallow water", long 
wavelength approximation in these derivations. 
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IV. IRROTATIONAL BACKGROUND FLOW 

We can use the equations of motion of the fluid to de- 
rive the most general rotationally symmetric and locally 
irrotational backgound flow profile. Let us assume that 
the bottom of the tank is defined by the relation 



Z = f(R) , 



(41) 



in the usual cylindrical coordinates (Z,R,ip), where / 
denotes some moderately curved function. Switching to 
the adapted coordinates (z, r, if) described in the previ- 
ous Section (Z = f(R) «-» z = ; r = R) gives the spatial 
metric 



dr 2 



dz 2 +dr 2 (l + /'(r) 2 ) 



f"(r) 



V(i + /'M 2 ) 3 , 



+ dip 2 



f'(r) 



for which the lowest order metric is 

dz 2 + (l + f'{rf) dr 2 + r 2 dip 2 . 



(42) 



(43) 



To lowest order in h the flow equations are the effective 
continuity equation (p0[) 



(44) 



dr(Vl + r(r) 2 rhv r ) =0. 



the condition for a locally irrotational flow 



(45) 



with L being some constant related to the angular mo- 
mentum, and finally the Bernoulli equation 



(l + /'(r) 2 ) K) 2 +r 2 K) 2 

h hoc cf s 

— Q — = — Qj( r ) , 

VT+7W 



(46) 



where we have neglected the "centrifugal" term in Eq. 
(E53|) as it will be very small for our situation. These give 



Ch Q 



rhy/l + f'(ry 



and 



1 C 2 h 



2 h 2 



L 2 

— I = 



h — h r 



V 1 + f'(rY 



(47) 



- 9f(r) . (48) 



Thus, we either need h to change as a function of r, or we 
need a non-trivial f(r). Choosing f(r) = —F/r 2 allows 
a consistant solution with contant height h — hoc for the 
fluid as long as F is given by F = (C 2 + L 2 )/g. The 
effective metric for the fluid is then of the form 



ds^ff = ^ ( ghao 5 ) dt z 



+2^ (yTTfWjdtdr + Ldtdipj 

-V([l + /'M 2 ]dr 2 +r 2 V) , 
9 



(49) 



with g = g z = g/y/l + f'(r) 2 denoting the effective grav- 
itational acceleration. 

In summary the analogy to a curved space-time and 
the concept of an effective metric can still be applied in 
the case of non-horizontal flow provided that the local 
variation of the height of the fluid and the slope of its 
bottom are sufficiently small. Nevertheless, the global 
changes may well be significant. 

A variation of h and g Zl the component of the force 
perpendicular to the bottom, does in general also entail 
a change of the local velocity \/g z h of the gravity waves. 
Such a spatial dependence may lead to further interest- 
ing effects: In analogy to optics one may introduce an 
effective index of refraction which then also acquires a 
non-negligible gradient. In such a situation the gravity 
waves may be scattered by this gradient or even the phe- 
nomenon of total reflection could occur. As we shall see 
later in Section [x|, this mechanism may be one ingredient 
for generating an instability. 



V. SURFACE TENSION 

So far we have considered ideal fluids without any in- 
ternal forces. However, if we take the surface tension 
of the liquid into account, the pressure at its surface no 
longer vanishes. Accordingly, the upper boundary condi- 
tion (||) for the pressure is modified to 



p(z = h) = — aV|/i , 



(50) 



where a denotes the fluid's surface-tension coefficient and 
V11/1 is the curvature of its surface in the linear approxi- 
mation. Consequently we obtain 



Spfts) — gg Sh — qV»(5/i , 



(51) 



instead of Eq. (|9[). This results in a extra term in the 
velocity perturbation equation (0) 



SQ + vl ■ V||<5* = -gSh+^-VjSh. 



(52) 



As we shall see below, the effects of surface tension be- 
come relevant for small wavelengths only. In this limit 
we may neglect the variation of the background flow 

V11 ® i>g w and obtain a modified wave equation 



d_ 

m 



5$ = .g/i B Vjj<5<I> 



Vjj<5$. (53) 
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In terms of the capillary constant a 2 = a/(gg) and the 
velocity of the unperturbed gravity waves c B = gh-e, this 
wave equation results in the following dispersion relation 



,11 



k 



2lA\ 



a z k 



(54) 



Therefore the incorporation of the effects of surface ten- 
sion leads to a "superluminal" dispersion relation (in the 

terminology of ]2^]) since (for v B = 0) the group velocity 
duj/dk as well as the phase velocity uj/k exceed cb for 
large wavenumbers k. 

However, we should bear in mind that the above calcu- 
lations are still based on the assumption of long gravity 
waves A ^> /ib- For short gravity waves A <C /ib> on the 
other hand, the dispersion relation reads (for = 0, cf. 

{§) 



gk(l 



,2,„2\ 



(55) 



Hence we can use the ratio a/h^ in order to alter the dis- 
persion relation for large wavenumbers k: For o»/ib the 
capillary waves dominate before the wavelength becomes 
smaller than the height, and we have a superluminal dis- 
persion relation, whereas for a <C Jib the short gravity 
waves dominate before the surface tension becomes im- 
portant, and thus one initially has a subluminal disper- 
sion relation, before the capillary waves finally take over 
at very short wavelengths. 

For example, for mercury the surface tension coeffi- 
cient a is about a m 0.46 N m _1 at room temperature 
293 K and hence its capillary constant a sa 1.9 mm. For 
water at 293 K we have a w 0.0725 N m -1 and hence 
a ~ 2.7 mm. This quantity can easily be manipulated 
by changing the temperature, adding surfactants, or by 
changing the fluid used. 



VI. VISCOSITY 

The dynamics of a viscous but still incompressible fluid 
are governed by the Navier-Stokes equations 

— =v + (v-V)v = -+g + vV 2 v, (56) 

at g 

where qv denotes the dynamic viscosity of the liquid and 
v its kinematic viscosity. 

The boundary conditions have to be modified as well. 
Instead of Eq. (^|) we have now 

p(z = h) = 2gvd z v z (z = h) , (57) 

and there are two additional restrictions on v\\ 

d a v\\(z = h) = -V\\v,(z = h), (58) 

and 

«u(z = 0)=0. (59) 



Let us investigate the effects of a finite but small vis- 
cosity on the wave propagation - where we restrict our 
examination to the case of a vanishing background flow 
i>b = for simplicity and employ the plane-wave ansatz 
with a frequency to and a wavenumber k. 

For an incompressible fluid the divergence of the lin- 
earized Eq. (Effl yields (for vb = 0) 



V 2 Sp 



(d 2 -k 2 )5p = ) 



(60) 



which has the solution Sp = Acosh(kz) + Bsmh(kz) with 
k = \k\. 

In the long-wavelength limit fcft-B <C 1 we may approx- 
imate sinh(fcz) w kz and cosh(fcz) »1+ (kz) 2 /2. Insert- 
ing the resulting expression back into the Navier-Stokes 
equations (B6h yields 



k — i— ) Sv z 
v , 



Ak 2 z + Bk 
gv 



(61) 



Defining k via k 2 — k 2 + iu>/v and SR(fc) > the general 
solution of this equation can be written as 



dVz = ^ Ak2 z + B k + Ce + kz + De -~ kz 
gvk 2 



(62) 



In addition to the long-wavelength limit A > Ab we 
assume the viscosity v to be very small v <C u>h B . In this 
case the exponentials exp(±fcz) are very rapidly varying 
functions - which can be used to simplify the analysis. 

Combining the equation of continuity d z 8v z +ik-5v^ = 
with the velocity boundary condition in Eq. (|5^) one 
obtains (d 2 + k 2 )5v z (z = Hb) = 0. In view of \k\ k 
and kfiB 3> 1 this implies that C is extremely small C oc 
exp(— kh,B)/k 2 and thus can be neglected. 

On the other hand, from 5v z (z = 0) = we obtain 
D = Bk/(gvk 2 ). Therefore, the term Dexp(— kz) in Eq. 
(p3) is relevant in a very thin boundary layer of order 
\J vjuj over the bottom only, cf. [M . 

In analogy to Eq. (^|) we may linearize the boundary 
condition for the pressure in Eq. (|5^) which fixes the 
integration constant A. 

The remaining condition vu(z = 0) = 0, i.e. d z Sv z (z — 
0) = 0, can be used to eliminate D and hence B. As 
one might expect, the solution for Sv z in the presence 
of a small viscosity displays only slight deviations (B 
is of order ^Jv) from the linear profile 8v z oc z used in 
the previous Sections - as long as one is well above the 
aforementioned boundary layer. 

Finally, Eq. (^), i.e. Sv z (z = hs) = icoSh, enables us to 
derive the dispersion relation 



gliBk' 



'- + %)• 



(63) 



Here one can read off the characteristic damping time r 
after which the viscosity effects become significant 
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(64) 



One observes that high frequencies are damped faster. 
This tendency becomes much stronger in the regime of 
short gravity waves where 



1 

vk 2 
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(65) 



holds, see e.g. p6[ . As a result, the at a first glance 
undesirable effects of viscosity can be utilized to damp 
out potential high-frequency noise and so single out the 
interesting (medium-wavenumber) instabilities by tuning 
v. 

For example, water at room temperature has a kine- 
matic viscosity of v w 1 mm 2 s _1 . Assuming a height 
/ib = 10 cm and a frequency uj = 1 Hz we infer from Eq. 
( |64| ) a damping time of the order of one minute - which 
still seems to be sufficient. 

The kinematic effect of the friction within the liquid 
could be diminished by selecting a fluid with a relatively 
high density g (e.g. mercury[]) in comparison with its 
inner viscosity. For mercury at room temperature v is 
given by v ps 0.12 mm 2 s . 

However, there is also another problem induced by the 
finite viscosity: the ansatz for the z-independent back- 
ground flow profile used in the previous Sections is not 
appropriate anymore since the fluid sticks to the bottom 
of the basin, cf. Eq. @. 

One way to solve this problem is to move the bottom 
so that its velocity is the same as that of the fluid 
at least in the interesting region, for example near the 
horizon. (One might also imagine manipulating the fluid 
near the bottom and thereby effectively simulating this 
motion.) After transforming into its rest frame the above 
calculation demonstrates that the solutions derived in the 
previous Sections are still a very good approximation. 

Without moving the bottom the boundary condition in 
Eq. (|59|) enforces a significant z-dependence of the back- 
ground flow profile. For instance a constant horizontal 
force (needed for maintaining the stationary flow) implies 
a parabolic flow profile d b oc z{2H-b — z). Unfortunately, 
it is not possible to cast the full wave-equation into an as 
tractable form as in Eq. (^) allowing for the identifica- 
tion of an effective metric in this situation (this problem 
is currently under investigation). The main obstacle is 
that the flow is longer irrotational. 

Nevertheless, for other scenarios one might be able 
to overcome this difficulty. If one injects a nearly z- 
independent stationary inflow (driven by a turbine, for 



t The use of a fluid like mercury has just another advantage 
since one would be able to detect the gravity waves very easily 
via considering the reflection at the fluid's surface in this case. 



example) on one side of the basin, then the flow will ba- 
sically remain irrotational throughout the basin - outside 
a thin Prandtl boundary layer at the bottom, cf. p6| . 

In view of the relatively large velocities involved and 
the supposedly small viscosity this layer may well be tur- 
bulent. However, by an appropriate preparation of the 
bottom's surface (e.g. dolphin skin effect) the indued drag 
can be diminished. 

Since the properties of the flow outside the thin bound- 
ary layer are nearly the same as in the case without vis- 
cosity and in view of the remarks after Eq. ( |62] ) one would 
expect that the basic properties of the gravity waves as 
discussed in the previous Sections are not drastically af- 
fected by a small internal friction in this case. 

Beside the scenarios described above there is also an- 
other - more exotic - solution for the viscosity problem 
conceivable: a superfluid does not stick to the boundary 
and its vorticity is quantized. But the necessity of using 
gravity waves in say liquid Helium makes the experimen- 
tal realization far more difficult than is desirable. 



VII. ENERGY 

There are two different concepts of a metric in our 
model: Firstly, the Minkowski metric determining the 
length- and time-scales within our laboratory; and, sec- 
ondly, the effective metric - which is experienced by the 
gravity waves only. 

These two concepts lead to two distinct notions of en- 
ergy: For time-independent external forces (and inviscid 
fluids) the Noether theorem demands the conservation of 
the total energy of the flow. In addition, assuming a sta- 
tionary background flow profile, we may find a conserved 
energy associated to the gravity waves: 

At the boundary dG of the basin, the normal com- 
ponent of the velocity v ■ n has to vanish. In terms of 
the velocity potential <E> this corresponds to Neumann 
boundary conditions 



v ■ n \ dG = 



n ■ V<f> \og= . 



(66) 



This enables us to accomplish a spatial integration by 
parts and in complete analogy to the 2+1 dimensional 
curved space-time one may derive a conserved energy 



E = / dE M T» v = 



(67) 



where £ M = d/dt denotes the Killing vector associated to 
the physical laboratory time and the spatial hyper- 
surface. For the effective metric in Eq. ( |16|) we obtain 
the energy density 



n — 77 



(68) 



of the perturbations 
also @). 



53? (o), i.e. gravity waves (see 
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We observe that this energy density contains negative 
parts inside the ergoregion \v^\ > y/ghn. This observa- 
tion points to the possibility of an instability associated 
to the ergoregion. (A positive definite conserved energy 
density would prove stability.) The total energy of the 
fluid including the background flow is of course always 
positive. 

Note that the energy conservation law derived above 
is violated if G has a hole, such as at the drain, etc. This 
problem, however, arises for real black holes as well. 



VIII. NON-ROTATING BLACK HOLE 

If we neglect the small slope f'(r) C 1 of the bottom 
(cf. the remarks at the end of Section III) in Eq. j49| ) 
the constant quantities g and can be absorbed by a 
simple rescaling and we arrive at 



d s lft = ( 9 h c 



(J2 , £2\ Q 

dt 2 + 2— dt dr + 2L dt dip 



-dr 2 



2 d V 2 



(69) 



If we take L — 0, we exactly recover a Painleve- 
Gullstrand-Lemaitre (PGL) type metric |2j| 

ds 2 cS = [c| - w 2 (r)] dt 2 + 2w(r) dt dr - dr 2 - r 2 dip 2 , (70) 



with c| = ghao and w(r) = ^/l + f'(r) 2 v r = C/r. As it 
is well-known, by means of the singular coordinate trans- 
formation 



dt — > dt = dt + dr ■ 



w{r) 



w 2 (r) 



(71) 



the stationary PGL metric can be cast into the static 
Schwarzschild form 



dsls = (4 ~w 2 )di' 



W 



T dr 2 - r 2 dip 2 . (72) 



Obviously the horizon occurs when w 2 — c| = g/iB, i-c 
when the velocity of the (radially) flowing fluid exceeds 
the speed of the (long) gravity waves V gh^. An inward 
flowing liquid w < simulates a black hole whereas an 
outward flow w > evidently corresponds to a white 
hole. The black hole branch can be used to observe the 
inversion of the centrifugal acceleration mentioned in 
the Introduction and, of course, the trapping of the waves 



IX. WHITE HOLE 

The white hole branch w > of Eq. (f70|) offers another 
interesting phenomenon: As demonstrated in Ref. all 
incident waves pile up at the horizon (since they cannot 
penetrate) and get arbitrarily blue-shifted there - if one 



neglects the change in dispersion relation, and thus group 
velocity, at high wavenumbers. In our model, however, 
the blue-shifted waves eventually leave the regime of the 
long gravity waves. The subsequent behavior depends 
on the character of the dispersion relation, in particular 
whether it is subluminal or supcrluminal. In the latter 
case the incident waves will eventually penetrate the hori- 
zon, once their wavelength has become sufficiently short 
to alter the group velocity to one larger than the (low 
frequency) velocity of the waves. In the former case, the 
short wavelengths will be swept out by the fluid which 
is now flowing faster at the horizon than the group ve- 
locity. These features associated with white holes may 
thus be observed experimentally for these gravity wave 
analogs. Furthermore, since the dispersion relation can 
be adjusted by varying the depth and the surface tension 
of the liquid, one can study the effect over a wide range 
of physical situations. 



X. ROTATING BLACK HOLE 

Stationary flow profiles containing a component in in- 
direction (vortex solutions with non-zero L) can be used 
to model rotating (Kerr) black holes. With a further 
rescaling and a redefinition of the constants C and L 
we may absorb the speed of the gravity waves ghoo in 
Eq. ( |69| ) completely and then the corresponding effective 
metric assumes the following form 



ds 2 eS = [ 1 



C 2 



L' z 



-dr 2 



r 2 d^ 2 



C 

dt 2 + 2—dtdr + 2Ldtd(p 



(73) 



The space-time structure of a Kerr black hole is more 
complicated than the Schwarzschild geometry for there is 
a difference between the static limit (or the ergosphere, 
see e.g. 0]) and the horizon: The static limit goo = 
denotes the region beyond which no particle can remain 
at rest. This, as mentioned above, corresponds to the 
surface where the velocity of the fluid equals the velocity 
of the waves, i.e. r 2 — C 2 + L 2 . The horizon is the 
"point of no return" , and for an axially symmetric flow, 
corresponds to the surface where the radial flow velocity 
equals the velocity of the waves r 2 — C 2 . 

The region between these two critical points, in both 
Kerr space-time and this model, is called the ergoregion 
and allows for the occurrence of the superradiant modes. 

According to Eq. ( |68| ) the energy density may become 
negative inside the ergoregio n r 2 < C 2 + L 2 . As al- 



ready anticipated in Section VTj , this observation can 
be interpreted as an indicator of instability. Indeed, in 
complete analogy with the Kerr black hole this analog 
should exhibits the phenomenon of superradiance: An in- 
cident wave with non- vanishing angular momentum scat- 
ters from the region around the black hole (analog) - i.e. 
the vortex - and the amplitude of the reflected wave is 
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larger than that of the ingoing wave. The necessary en- 
ergy is extracted from the rotational energy of the back- 
ground. 

Since the metric in Eq. ( |73| ) possesses two independent 
Killing vectors, d/dt and d/dip we may find a complete 
set of solutions of the wave (KFG) equation 







2- d t d r 
r 



2^ d t d v + — ^ d r d v 



-9 r — a, 



L 2 

r 4 



1 



d 2 



by the following separation ansatz 

(j>(t, r, <p) = exp {— iut + irmp} <p um (r) . 



= 0, (74) 



(75) 



Note that an analogous separation ansatz in t, i?, <p is even 
possible in the real 3+1 dimensional Kerr metric - which 
is a less trivial statement. 

The remaining function <f>w m (r) obeys a second-order 
ordinary differential equation. In terms of the Regge- 
Wheeler tortoise coordinate defined by 



the KFG equation (fj 
— oo, simplifies to 



dr . 



r 2_ C 2 

\) at the horizon r 2 



(76) 



C 2 , i.e., r* 



8 r 



2i 



C 2 



9 r 



0. 



(77) 



The two linearly independent solutions to this equation, 
i.e., (approximately) constant and purely oscillating, re- 
spectively, correspond to ingoing and outgoing modes, 
respectively. 

In analogy to the Kerr metric we may introduce the 
angular velocity of horizon 



L 



(78) 



For a Kerr black hole this quantity is bounded]]] by its 
total mass M via 20h < 1/M. With the gravity wave 
analogs, however, it is possible to generate rather large 
values of 2Qn - in particular if one allows for non- 
negligible slopes - since C and L can be varied inde- 
pendently. 

The Wronskian associated to the second-order ordi- 
nary differential equation reads at the horizon 

w wm [0] = 0*a P .0-0a P .0* + 2i(w-n H m)^. (79) 

Inserting the ingoing mode d r , (j) — Q we obtain 



* Otherwise the metric would describe a naked singularity 
without ergoregion and horizon, etc. 



2i(u> — flftm) 



(80) 



where |7^ m | 2 = |</> wm | 2 denotes the transmission coeffi- 
cient. Since we have restricted our solution to be purely 
ingoing at the horizon, it will contain ingoing as well 
as outgoing components (determined by the reflection 
coefficient lZu, m ) at spatial infinity r — oo in general 
<t> oc exp(— iu)r) + K ura exp(io;r). In this limit the wave 
equation (|74|) reduces to the usual form and hence the 
associated Wronskian reads 



W um [4>] = 2kj(1-\H 



(81) 



For a regular second-order ordinary differential equa- 
tion on a contiguous interval the Wronskian is conserved 
which implies Q the following relation between the trans- 
mission 7^ m and the reflection K ura coefficients 



■ |2 to — mfin ^ 12 



(82) 



As one can easily infer from the equation above, if 
mQn > u> the reflection coefficient is greater than one 
l^wml > 1) i- e - t ne scattered wave has a larger ampli- 
tude than the incident one. This amplification process 
corresponds to the phenomenon of superradiance. 

However, for a massless scalar field in the asymptot- 
ically flat Kerr geometry the scattered wave escapes to 
infinity. To make an amplifier unstable, feedback is re- 
quired, see p(i|| . That outgoing wave must be reflected 
back toward the hole for repeated amplification. This is 
possible if the field has a finite mass |J or if one, for exam- 
ple, encloses the rotating black hole by a large spherical 
mirror 101. 



As explained in Section IV, for the gravity wave 



analogs there is another mechanism which may force the 
scattered wave to "come back" - the total internal reflec- 
tion. Since the velocity of the wave is \gh, if one can 
manipulate h, as happens for example near the drain of a 
bath tub, the change in the effective refractive index may 
entail effectively trapping the waves - which then can be 
amplified via superradiant scattering without bound, or, 
rather, until non-linear effects dominate. 

In summary it is possible that some of the features 
which one regularly observes in the vortex flow near the 
drain of one's bathtub are in fact the result of instabilities 
which are exactly analogous to the behaviour of waves 
near a (rotating) black hole. It should be mentioned here 
that the arguments above do not explain the formation 
of the vortex, but refer to its linear instability. 



XI. INNER HORIZON 

The instability indicated by negative parts of the en- 
ergy density in Eq. (|6^) is not necessarily restricted to 
rotating black holes: In the presence of an inner hori- 
zon in addition to the outer one, such as occurring in the 
Reissner-Nordstrom metric, and for an altered dispersion 
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relation, runaway solutions can exist even in the purely 
one dimensional flow. 

This somewhat surprising fact has been demonstrated 
in Ref. ^7j. In the following we give a brief repetition 
of the basic explanation adapted to subluminal high fre- 
quency dispersion relation. Note that the terms "sub-" 
or " superluminal" do not, despite their origin, refer to 
light, but rather to whether or not the group velocity at 
high wavenumbers is less than or greater than it is at 
very low wavenumber (long wavelength). They will also 
refer to regions where the velocity of the fluid is smaller 
or larger than than the low-wavenumber velocity of the 
waves. 

We shall restrict ourselves to a one-dimensional flow. 
We consider a case where there are two horizons, one a 
white hole (where the velocity of the fluid drops below 
the low wavenumber velocity of the wave) and the other a 
black hole (where the velocity of the fluid goes above the 
low wavenumber group velocity) . It turns out (for not en- 
tirely well understood reasonsjj) that although there are 
two branches of the dispersion curve (for a static fluid 
these correspond to the two directions of propagation of 
the wave), wave packets remain on one or the other of 
these branches even when interacting with horizons, and 
spatially varying flows. There seems to be very little 
backscatter from one branch to the other. We shall re- 
strict attention to the one branch which at low wavenum- 
bers represents waves which are travelling in a direction 
opposite to the fluid flow. 

Assuming adiabatic motion of this wave packet (in 
analogy to the geometric optics approximation), the ef- 
fect of the fluid flow is to alter the dispersion relation 
to 



duo 
dk 



dF 
~dk' 



(84) 



which is negative, corresponding to travel to the left, 
against the flow of the current. The other two, one at pos- 
itive fl and one at negative Q both have positive group 
velocity and thus correspond to packet travel with the 
current flow. (Their rest frame group velocity is so small 
that even the subluminal current is sufficiently fast to 
drag them along) . 
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to + vk = F{k) , (83) 

where the mode has the form e l ( UJt + kx ) ^ i. e . 5 u i s the 
frequency of the wave in the lab frame, not the rest frame 
of the fluid. Define Q = lo + vk as the frequency in 
the fluid rest frame. In Fig. || we sketch the subluminal 
dispersion relation as a plot of f2 versus k. To determine 
the possible values of k for any given value of lo, we can 
plot the line Q — vk = lo and look for the intercepts 
with the dispersion curve. We shall concentrate on small 
values of lo, which is also the intercept of the line with 
the ft axis. 

On the subluminal side of the horizon, where v and 
thus the slope of the line is less than the dispersion slope 
at k — 0, there will in general be three points of intersec- 
tion, and with very different group velocities. The one 
near the origin k = has a group velocity, 



^The 1+1 dimensional scalar wave equation without disper- 
sion is conformally invariant which explains the absence of any 
mixing in this case. For a non-trivial dispersion, however, this 
problem is less clear. 



FIG. 2. One branch of the dispersion relation where the 
high wavenumber group velocity is smaller than the low fre- 
quency. Plotted are the two lines of Q vs. k for different 
velocities, but the same value of the lab frame frequency lo. 
Note the intersection points of these straight lines with the 
dispersion curve indicating the possible values or k for the 
given lo at the given fluid velocity. 

Because the background current flow is assumed to be 
stationary, the frequency lo is conserved during the mo- 
tion of a wave packet. Let us consider what the possible 
wavenumbers are for this same value of lo on the super- 
luminal side. The slope of the line is now much greater 
than the low wavenumber slope of the dispersion curve, 
and there is now only one intersection with the dispersion 
curve, at negative values of f2. Furthermore the group ve- 
locity here is positive. These wave packets, even though 
they have low wavenumber, are dragged along with the 
fluid and travel in the same direction as the fluid. 

One can, as for scalar waves, define a "norm" of 
the various wave packets by the conserved Klein-Fock- 
Gordon inner product 
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($1$) =i J d^(<P*d^~$d^*) 

= l -Jdx<S>* (dt -v{x) d x ) $ 

= Cl(x) J dx\<$>\ 2 , (85) 

where we have assumed the wave packet to be sufficiently 
localized that v is constant over the packet. I.e., the sign 
of the norm is the same as the sign of f2. Thus, for 
positive uj the wave packets on the superluminal side will 
have negative norm, while those on the subluminal side 
can be either positive or negative norm (two are positive 
and one is negative). 

Let us consider a wave paket with positive norm, low 
value of uj but large value of k being dragged toward 
the black hole horizon. As it comes closer, v increases, 
and the wavelength of the packet is stretched out, with 
k becoming smaller and smaller. As it hits the horizon, 
where the slope of the line becomes essentially the same 
as that of the low wavenumber dispersion curve, the adia- 
batic assumption (geometric optics approximation) fails, 
and one gets a mixing of the various possible values of k 
for the given u. However, whatever happens, the wave 
packet must thereafter leave the horizon, either on the 
subluminal side, or the superluminal. On the subluminal 
side the only possibility is that the wave exit with the 
small value of k near zero, as this is the only one of the 
three possibilities which travels away from the horizon on 
the subluminal side. On the superluminal side, the only 
possibility is that the wave packet leave with negative 
norm. The fluid flow at the horizon will have mixed pos- 
itive with negative norm solutions. This is precisely the 
requirement in the quantum system that particle produc- 
tion take place, since this mixing of positive and negative 
norms is precisely what leads to non-trivial Bogoliubov 
coefficients 



negative norm packet into the superluminal region, al- 
though if the phases were just right, the two could cancel 
and emit nothing into the superluminal side. Assuming 
that this does not happen, or that the conversion at the 
white hole is so small so as not to create any mixture, 
the subluminal packet's norm increases once again. Af- 
ter many such back and forth reflections, the norm can 
grow arbitrarily large. One has an instability. In Fig. ^ 
we have four steps in this process. The black hole horizon 
occurs in the center of the diagram, while the white hole 
horizon is at the left and right edges (we assume periodic 
coordinates). The low frequency intial wave packet trav- 
ells toward the white hole horizon, and reflects as a high 
frequency wave, dragged by the fluid. The third window 
is after the reflection from the back hole, and the fourth 
after the refection again from the white hole. 
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P = (#LI*out) • (86) 

For our purpose, what this means is that the wave packet 
which leaves the horizon on the subluminal side must 
have a larger norm than did the wave packet which was 
dragged toward the horizon. (Since the total norm is 
conserved, and since the norm of the packet on the su- 
perluminal side is negative, the subluminal packet must 
be larger.) 

This packet will now travel toward the white hole hori- 
zon. Here, there is no packet on the superluminal side 
which travels away from the horizon. The only possibility 
is that the packet be blue shifted at the white hole horizon 
and come off as a mixture of the positive and negative 
large k solutions. While one of these is positive norm, 
and the other is negative, they are travelling together 
toward the black hole horizon. Here they again emit a 



FIG. 3. Four stages in the process of a wave packet re- 
flecting off the white and black hole horizons. The black hole 
horizon is in the centre. These are the real parts of the wave 
packets (the imaginary look similar and were chosen to make 
the intial pulse purely positive norm). 

For small oj the amount of mixing of positive and neg- 
ative norm modes at each of the horizons is governed by 
a thermal Bose-Einstein factor 

l/3| 2 « 1 , (87) 

exp(o;/i ) — 1 

with an effective temperature T proportional to the effec- 
tive surface gravity, i.e., the rate of change of the velocity 
or fluid flow across the horizon. If one makes the white 
hole horizon have a very low effective temperature, one 
can minimize the creation of large k negative norm so- 
lutions on "reflection" of the packet from this horizon. 
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This means that the packet on the subluminal side will 
remain positive norm, if it started thus. 
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FIG. 4. The reflection of a pulse producing low frequency 
instability. Plotted are three steps in the evolution of a pulse 
(top box) after one (middle box) and two reflections from the 
black hole horizon in the center. Note that in each case the 
pulse on the left (in the subluminal region) is travelling to the 
left while the ones on the right (in the superluminal region) 
travel to the right. After another reflection from the black 
hole horizon, the pulse grows linearly and after still another, 
quadratically with distance from the black hole horizon. 

One of the consequences of this instability is that one 
expects it to be worst for very small values of u>. Since 
the norm mixing at the black hole horizon is proportional 
to a Bose-Einstein factor, this diverges for small values 
of ui 



exp(w/T) - 1 to ' y J 

One would thus expect, and numerical simulations con- 
firm, that after a few back and forth bounces off the 
black hole horizon, the lowest wavenumber mode of the 
wave (i.e., constant in space) would be amplified at each 
bounce. Numerical simulation shows that after a couple 
of bounces, at each reflection from the black hole horizon, 
a constant "step function like" pulse of ever increasing 
amplitude is emitted by the black hole horizon into the 
subluminal side, cf. Fig. |[ After the next reflection the 
spatial groth is linear, then quadratic etc. 

One can, following Corley and Jacobson p7j , make a 
very similar analysis in the case in which the dispersion 
relation is superluminal 



F(k) = c 2 Vk 2 + a 2 k i , (89) 

which is exactly the situation considered in Section In 
that case one should consider a superluminal region be- 
tween a black hole and a white hole horizon. In this case 
it is the wave packet bouncing between these two hori- 
zons on the superluminal side which creates the instabil- 
ity. Taking positive co modes but with an intial negative 
norm wavepacket, we find that every time the packet hits 
the black hole horizon it emits a positive norm packet 
into the subluminal region beyond that horizon (which is 
again exactly the classical counterpart of Hawking radia- 
tion), thus increasing the size of the negative norm packet 
in the superluminal region, while every time it bounces 
off the white hole horizon, it simply bounces off, possibly 
mixing in some high frequency positive norm packet as 
well (depending again on the effective temperature of the 
horizon). 

Since with gravity waves, we can choose whether to 
have a subluminal or superluminal dispersion relation, by 
choosing the depth and surface tension of the liquid, one 
can test both sets of predictions. Furthermore, since this 
norm mixing property of the horizon is intimately related 
to the quantum emission of radiation by the horizon, one 
can, with purely classical gravity waves, test the predic- 
tions about the thermal behaviour of the various types 
of horizon. 

The transition between short and long gravity waves 
can be described by the height-dependent dispersion re- 
lation for a fluid at rest @ 

lj 2 = gkt&nh(kh B ) = gh B k 2 ^1 - + 0(k 6 ) . (90) 

This should be contrasted with the dispersion of the sur- 
face tension waves as derived in Eq. ( |54| ) of Section |v| 
(again for a fluid at rest) 

lj 2 = g h B k 2 (l + a 2 k 2 ) . 

Consequently, one obtains a superluminal dispersion re- 
lation for a ss h B which implies h B ss 2 mm for mercury. 
For a <C h Bl on the other hand, the dispersion relation 
possesses a large subluminal region (short gravity waves) 
before it becomes superluminal (surface tension waves) 
again. 

A possible experimental set-up for observing this in- 
stability is sketched in Fig. The fluid is flowing within 
a circular basin with varying depth such that the liquid's 
flow velocity as well as the speed of the gravity waves 
depends on the position. As the depth of the basin de- 
creases the velocity of the fluid increases and in the same 
time the speed of the gravity waves decreases. 

In this way one may construct a region of superluminal 
fluid velocity confined between the inner and the outer 
horizon. As mentioned above, the gradient of the basin's 
depth (determining the the surface gravity) should be 
small at the inner (white hole) horizon whereas it should 
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be as large as possible - without violating the assump- 
tion in Section e.g. by generating a breakdown of the 
laminar flow - at the outer (black hole) horizon in order 
to observe the instability described above. 

On the other hand, one also has to make sure that the 
wave packets emitted by mode conversion into the other 
region do not travel around and disturb the system when 
they come back through the other horizon. For real black 
holes they escape to infinity. In our system, however, this 
is not the case and thus one has to install an additional 
absorbing device, which effectively models spatial infinity 

in order to prevent this interference. 

superluminal 
flow 



black-hole 
horizon 



white-hole 
horizon 




absorbtion device 
(spatial infinity) 

FIG. 5. Sketch of a possible experimental set-up. A 
wavepacket is bouncing back and forth in the superluminal 
region between the inner (white hole) and the outer (black 
hole) horizon, emitting a piece of "Hawking radiation" at the 
black hole horizon and thereby increasing its amplitude. 

According to the remarks in Section [Vl| it appears to be 
reasonable to move the bottom (e.g. made of rubber) of 
the fluid or exert a magnetic force on it. In practice one 
might generate a flow with a slowly increasing velocity 
until the two horizons are formed. In this case one should 
observe an instability of the laminar and smooth flow 
profile exactly at this threshold velocity. Of course one 
has to ensure that there is no other instability at or before 
this thershold which might spoil the experiment. 



XII. DISCUSSION 

In view of the main advantage of the gravity wave 
analogs - the possibility of tuning the velocity of the 
wave propagation rather independently - one would ex- 
pect that it is possible to simulate black and white holes 
in the laboratory and to study their instabilities exper- 
imentally. As one perhaps reasonable hierarchy of the 
different dimensions involved 5h <C *C A one could 
imagine waves of about 1 mm amplitude in a 1 cm deep 
basin with the wavelength being circa 10 cm and the 
characteristic size of the black or white hole analog (its 
Schwarzschild radius) approximately 1 m. The velocity 



of those waves - and hence also that of the fluid - in this 
case would be about 0.3 ms" 1 and should be realizable. 
Aiming for the incorporation of surface tension waves 
for large wavenumbers it might be suitable to divide the 
above suggested length scales by a factor of about 10. 
Here we encounter another advantage of the gravity wave 
analogs: In contrast to sound waves, for example, the am- 
plitude of the gravity waves can be measured directly (as 
a length) and with a very high accuracy (e.g. via inter- 
ferometry using the reflection at the fluid's surface). 

However, all the suggestions above have to be recon- 
sidered when aiming for an experimental verification of 
the Hawking effect |ll| . In order to detect this quantum 
instability the temperature of the fluid should be as small 
as possible - whereas the fluid's velocity and its gradi- 
ent have to assume their maximal feasible values. The 
experimental verification of the Hawking effect seems to 
be on or even still beyond |2j| the edge of our present 
experimental capabilities - for the dielectric/optical as 
well as for the sonic/acoustic black hole analogs (Dumb 
holes) and probably even more so for the gravity wave 
analogs. 
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